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The lasing threshold is studied of a multi-mode chaotic cavity (linear size D 3> wavelength A) 
coupled to the outside through a small hole (linear size d <C A). For sufficiently weak absorption by 
the boundaries, the statistical distribution of the threshold is wide, its mean value being much less 
than the pumping rate needed to compensate the average loss. The average number (Nnc) 1 of 
non-competing excited modes is proportional to the square root of the pumping rate. We use the 
classical model of spatial hole burning to account for mode competition and find a reduction in the 
average number of excited modes to (N) = 

PACS numbers; 42.55.Ah, 05.45.-fb, 42.55.Sa, 78.45.-bh 


I. INTRODUCTION 


Incorporation of quantum optical effects is a necessary 
and interesting extension of active ongoing research on 
multiple scattering of electromagnetic waves in random 
media ||^. It becomes particularly important when the 
medium is active, as is the case in experimentally realized 
“random lasers” Quantum effects have been largely 

ignored in many publications devoted to propagation in 
disordered amplifying waveguides in which only 

amplified stimulated emission of external incoming flux 
but not of spontaneously emitted internal noise was taken 
into account (Ref. being a notable exception). Ampli¬ 
fied internal noise can lead to excitation of low-threshold 
lasing modes of the waveguide, making practical use of 
amplifying waveguides problematic. The difficulty of the 
waveguide geometry is the onset of localization. In this 
paper we consider a simpler cavity geometry, which does 
not show localization, but retains two essential features 
of the problem: large sample-to-sample fluctuations and 
spontaneous emission. 

A complete description of the fluctuations is possible 
in the universal regime characterized by a chaotic pattern 
of classical trajectories. We assume that the cavity (vol¬ 
ume V ~ D^) is confined by conducting walls, filled with 
a lasing medium (central frequency of the gain profile 
Wo), and coupled to external detectors via one or several 
small holes. It was demonstrated recently that a non- 
integrable shape of the resonator can significantly affect 
its lasing properties Chaoticity of classical trajec¬ 
tories can be achieved either by a peculiar shape of the 
resonator Ml , or by a small amount of disorder scat¬ 
tering. We will speak about “chaotic cavities,” meaning 
any of the two mechanisms responsible for the onset of 
chaos. 

We restrict ourselves to the case of well-resolved cav¬ 
ity modes, which means that 1) resistive loss 7* in 
the cavity walls is less than the mean modal spacing 


( 5 wo = tt^c^/wqU, 2) characteristic size of the holes d 
is smaller than the wavelength Aq = 2 ttcIujq. Mean loss 
7o through a small hole was calculated by Bethe |M|, 


70-^^, d<Ao, (1) 

so that 7 o/(5wo — (d/Ao)® 1 . Note that the loss (||) 

is not proportional to the area of the hole. It is in fact 
much smaller than one might guess by extrapolating the 
dependence 70 — cd‘^/V valid for d ^ Aq. The effect 
of sample-to-sample fluctuations is pronounced only if 
7* 7o. This regime is experimentally accessible, as 

was demonstrated by a recent series of experiments on 
microwave cavities with superconducting Niobium walls 



Each act of spontaneous emission in a pumped cavity 
is a source of radiation into some cavity mode. Classical 
condition of the lasing threshold in a given cavity mode is 
satisfied if the gain due to stimulated emission equals the 
loss. Threshold for the cavity is the smallest value of the 
pumping rate at which threshold is attained for one of the 
modes. The questions we ask are: What is the threshold 
rate of pumping? How many lasing modes can coexist for 
a given pumping rate above the threshold? The problem 
of spectral content of outgoing radiation has been widely 
studied for integrable cavities of definite shape. Consid¬ 
ering arrays of chaotic cavities of slightly varying shape 
or with different configurations of scatterers we address 
the problem statistically and compute the probability of 
lasing, the distribution of the threshold, and the average 
number of excited modes. 

Trivially, gain greater than mean loss 70 will be on 
the average sufficient to ensure lasing, while gain smaller 
than 7* will never suffice. The mean loss from a tiny 
hole is small. We argue that the actual average thresh¬ 
old can still be many orders of magnitude smaller. Each 
individual cavity exhibits a well-defined threshold but its 
statistical distribution is wide. In Section Q we compute 
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this distribution for the idealized case 7 * = 0. Effects 
of non-zero resistivity of the walls are discussed in Sec¬ 
tion III. Section IV is devoted to the computation of the 
average number of excited modes above the threshold. 
We conclude in Section 


II. DISTRIBUTION OF LASING THRESHOLD 

We assume that the line of spontaneous emission is ho¬ 
mogeneously broadened and has Lorentzian shape with 
central frequency wq and width 20. Let ipi{^ be the am¬ 
plitude of a mode of the closed cavity at frequency LOi, 
normalized according to jdriplif) = V. (For simplic¬ 
ity we neglect polarization dependent phenomena and 
work with real scalar field amplitudes.) In the presence 
of weak coupling to the outside world the modes acquire 
finite widths, 7 ^. We assume two sets of conditions: 

70 < (5wo < O < Wo, d < Ao < U. (2) 

An especially important role is played by the inequality 
7 o SojQ, which is implied by d <C Aq. It ensures that 
the modes of the open cavity are well-defined and do not 
differ significantly from those of the closed one. In this 
Section we consider the idealized case in which there is 
no loss in the walls of the cavity ( 7 * = 0 ). 

In a chaotic cavity the modes can be modeled as 
random superpositions of plane waves [^. (Validity of 
this model has been checked experimentally in microwave 
cavities This implies a Gaussian distribution 

for ipiif) at any point r. The corresponding distribu¬ 
tion for '!/'?(r) is called the Porter-Thomas distribution 
[ p^ . Loss from a small hole located at f is proportional 
to [Vs' 0 i(r)]^ (with the derivative in the direc¬ 

tion normal to the surface of the hole) and has the same 
Porter-Thomas distribution, which was directly probed 
in the experiments of Ref. |]^. More generally, the dis¬ 
tribution of normalized modal widths y = 7^/70 in a 
cavity with v holes is given by the y^-distribution with 
v degrees of freedom (normalized to 1 ), 

Pu{y)= p(^/2) ?/~^^''^^exp(-i^y/2). (3) 

We assumed that loss from different holes is indepen¬ 
dent, which is true provided their separation is larger 
than Aq. For small integer v, the distribution (^) is wide. 
The single-hole case v = 1 looks especially promising 
from the point of view of low-threshold lasing because 
Piiy) = exp(—j// 2 )/-y/ 2 ^ grows with decreasing y. 

To grasp the picture we first confine ourselves to a 
subset of cavity modes located near wq. We neglect fluc¬ 
tuations of their frequencies and assume that the modes 
are equidistant, Wm = wq -I- Stuom, m = 0,±1,±2,... 
We denote by Rpg a reference pumping rate necessary to 
provide gain equal to the mean loss 70 at frequency wq. 


and introduce the reduced pumping rate e = Rp/Rpg, 
assumed 1. Loss of different modes is uncorrelated 
and distributed according to Eq. (^) while gain dimin¬ 
ishes with increasing difference |w — WqI according to the 
Lorentzian 

5o(w) =7oe[l + (w-wo)^/f2^] ^ (4) 

It follows that the probability Piy{e) of there being no 
lasing mode at the pumping rate e is given by 

ga{^m)/lo \ 

dyPAy)]- (5) 


pas ) =n (^ “ 


For e <C 1, the upper limit of the integral is also <C 1, 
and we can replace Pu[y) by its leading behavior at small 
y, Pv{y) oc which yields 

« exp X/ , ( 6 ) 

a = (z./2)-i+"/2[p(^/2)]-i. (7) 


Because the summand decays as we find that for 
V > \ the leading behavior of the probability of no-lasing 
is determined by the modes with |to| < Vt/dojQ, 


Pu{e) 


exp 


-C,{n/5LOo)e’^/^ 


( 8 ) 


Here Cl, = yCTC,yT{{h'—l)/2)/T{v/2) foT v > 1 (below we 
will separately define Ci). Modes far from wq have neg¬ 
ligible chance to get excited and need not be taken into 
account. On the contrary, for ^ = 1 all cavity modes, 
including those very far from wq, contribute to the prob¬ 
ability. 

To treat the contribution of distant modes for u = \ 
correctly, we must account for several factors which we 
could ignore for jz > 1. 1) The spectral density cannot 
be replaced by its value po = I/^wq at oj = loq. Instead 
p(w) = juA- 2) The mean loss is frequency depen¬ 
dent, 7 ( 0 ;) = 7 oa;^/wQ, cf. Eq. (|^). 3) The Lorentzian (^) 
for the amplification rate is an approximation valid only 
in the vicinity of wq. A correct expression for the gain 
g(uP) must be even in lo to comply with the symmetry 
x(w) = x*(—w) of the dielectric susceptibility y. It in¬ 
cludes contributions of both poles iwo + and reads 




4u;^7o£U^ 

[up — Wg)^ -I- 2[uj^ -I- Wg)!!^ -I- 


(9) 


Taking these three factors into account and replacing 
the discrete sum by an integral, the probability of no¬ 
lasing is given by 
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p^{e) = exp 



div p{io) In 


/ dyP^(y) 


. ( 10 ) 


For ly > 1 this leads to Eq. (^, the ultraviolet cutoff 
uJmax being irrelevant. For z/ = 1 we get 


Pi (e) = exp 



•.X 

duj p{iyj) In 


1 — erf. 


9{^) 

27 (w) 


( 11 ) 


where erf( 2 ;) = ( 2 /-y/ 7 r) dxexp(—x^) is the error func¬ 
tion. The main logarithmic contribution of type f duj/u) 
to the integral in Eq. ( 0 ) comes from large values of w. 
The ultraviolet cutoff iOmax — ‘i'Kcjd appears because loss 
of high frequency modes with A < d no longer exhibits the 
strong fluctuations of Eq. i)- Beyond the cutoff classical 
ray optics applies, leading to a narrowly peaked distribu¬ 
tion of the loss around the value cd^/E » 79 . Because 
we are considering the case £ <C 1 in which the gain 
is much smaller than the average loss 70 , the high fre¬ 
quency modes cannot be excited. It follows that the only 
relevant cavity modes are those with frequencies smaller 
than ujmax- Their number M u^^^/iujQSoJo — (D/d)^ 
is ^1. From Eq. & the probability of no-lasing pi (e) 
can be cast in the form of Eq. (||) with the coefficient 
Cl = ln(a;maa;/f^) weakly dependent on the fre¬ 

quency cutoff • Fig. 0 shows that the probability of 
lasing 1 —pi(£) can be reasonably large even for extremely 
small values of the reduced pumping rate £. 

The quantity 1 — Pv[e) is the fraction of lasing cavi¬ 
ties in an array at a given pumping rate e. It is directly 
related to the probability distribution Tv{£) of the las¬ 
ing threshold. Obviously dPT^^e') = 1 —pv{£), hence 
Tv{£) = —dpj/(e)/d£. We find from Eq. (||) that 


Tu{e) 


\vCi,{^/5ujD)e 


l+v/2 


exp 


-a(r!/(5a;o)£"/2 



FIG. 1. Probability of lasing 1 — pi(e) versus reduced 
pumping rate e with pi{£) given by Eq. ( |ll[ ) {loq/Q, = 10, 
9l/SuJo = 10). Thick lines are for different ratios D/d corre¬ 
sponding to different numbers M ~ (D/d)^ of relevant cav¬ 
ity modes (dot-dashed line M = 10^, dashed line M = 10^, 
solid line M = 10“^). Left inset shows an example of chaotic 
cavity. Chaotic behavior of classical trajectories in this par¬ 
ticular “die” shaped cavity was shown in Ref. 0 . Radia¬ 
tion is confined inside by means of ideally conducting walls 
and can leave the cavity only through a tiny hole. Right 
inset shows the probability distribution of the lasing thresh¬ 
old T^{x) = (v/2)x~^^''C with x related to e by 

X = e{Cv9l/dijjof'C ^ for different number of holes = 1, 2, 3. 


III. EFFECTS OF NON-ZERO WALL 
RESISTIVITY 


A non-zero loss 7 * from the resistivity of the cavity 
walls modifies the functions ( 00 ) by suppressing las¬ 
ing for £ < 7 */ 7 o- The distribution of the lasing thresh¬ 
old remains wide, as long as 7 */ 7 o “C I, as we now show. 
Instead of Eq. (0 we have 




pj.(£)=expl / dio p{Lo)\n j dy P„{y) 

\Juj- J(g(aj)-7.)/7(aj) , 


(13) 


( 12 ) 


(Deviations which arise at £ ^ I are unimportant.) 
The distribution is wide and in the single-hole case 
1 / = I diverges as £ ^ -bO (see right inset of Fig. |^). 
The average reduced threshold reads (e^) = r(I + 
2/z/)((7i,f2/^a;o)“^^'^. It is smallest for ly = 1 and is in¬ 
deed much smaller than I. 


where uj- < are the two positive frequencies such that 
g{uj±) = 7 *. A non-zero value of 7 * reduces the relevant 
frequency range to a narrow window around wq- There¬ 
fore, the modifications I)-3) of the previous Section be¬ 
come unnecessary even for the case = I. Using the sim¬ 
ple Lorentzian (^) for g{u!), instead of the more compli¬ 
cated expression (^, we find uj± = u;o±f^(£ 7 o/ 7 * — 1 )^^^- 
Neglecting the az-dependence of p{uj), 7 ( 0 ;) and using 
the small-argument behavior of the probability function 
Pu{y), we reduce Eq. (p^ to 
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p^{e) = exp\-C^{n/6ujo){-j^./joy^'^fui£lo/l* - 1 ) 


( 14 ) 


where Ci, is the numerical coefficient introduced in 
Eq. and 

can be expressed in terms of a hypergeometric function. 
In Fig. 1 ^ we have plotted the distribution of the lasing 
threshold, Ty{e) = —dp^{e)/de, for = 1 and different 
values of 7*/7o- We will analyze two limiting regimes. 



0 12 3 4 

EYod, 

FIG. 2. Probability distribution of the lasing threshold in 
a cavity with small absorption in the boundary and a sin¬ 
gle hole {u = 1), computed as —dpi{s)/ds from Eq. (^). 
We chose fl/Stoo = 10 and took three values of 7 */ 7 o such 
that 7*/7o is much smaller, equal, or much greater than 
{Sujo/O.f^'' = 10 "^ 


In the regime eyo /"f* S> 1 and for > 1 we recover the 
expression with the same constant C^. The value of 
Cl = Cl ln(e7o/7*) is different because of the different 
cutoff mechanism. Instead of having a weak logarithmic 
dependence on Umax it exhibits a weak logarithmic de¬ 
pendence on the pumping rate e. This limiting case is 
statistically dominant if 7*/7o , because 

then the corrections to Eq. at £ < 7*/7o have negli¬ 
gible statistical weight. 

In the opposite regime, £70/7* — 1 1 , the threshold 

distribution differs significantly from Eq. 


Ty[e) = i(l-bi^)A^(II/5wo)(7o/7H-)^^^(^-7*/7o)"=' 


X exp 


-Aj.(II/^a;o)(7o/7*)^^^(£ -7*/7o)'^" 


( 16 ) 


(with a numerical coefficient -|- 

i^/2)/r(3/2 + i^/2)). This regime is statistically domi¬ 


nant if 7*/7o ^ (< 5 a;o/f 2 )^/'^. The mean value of thresh¬ 
old is now close to 7*/7o, but there are large fluctuations 
towards larger e. 


IV. AVERAGE NUMBER OF EXCITED MODES 


In this Section we focus on the number of lasing modes 
beyond the lasing threshold for = 1 assuming 7 * = 0 . 
We assume that the parameters are such that many 
modes are above the threshold. This requires, in par¬ 
ticular, £ 70 / 7 * 3> 1. In this case a non-zero value of 7 * 
only leads to a redefinition of Ci because of the different 
cutoff mechanism. If the modes did not compete we could 
compute the average number of excited modes {Nnc) as 


{Nnc) 



du p{oj)eTi. 


g(‘^) 
2 j{u) ■ 


( 17 ) 


For £ < 1 it is given by (Nnc) = C'i(n/ 5 u;o)£^/^. How¬ 
ever, the modes do compete for a homogeneously broad¬ 
ened line because one of the modes can deplete the in¬ 
version, preventing another mode from being excited [ p^ . 
Multi-mode operation is still possible if different excited 
modes deplete the inversion in different spatial regions of 
the cavity ]l 9 | , p 0 | . We assume this mechanism of multi- 
mode generation, called spatial hole burning . 

Let ni, Af {r) denote the number of photons in the mode 
i and the density of population inversion between the las¬ 
ing levels. Semiclassical rate equations read 

^ = -7i’^» + W.rii Jdf' 4 ’‘f{r)Af{r), ( 18 ) 

= sRpolV - u)Af{r) - Af{f) ^ WiTiii^iir). ( 19 ) 


Here w is the non-radiative decay rate and is the 
rate of stimulated emission into mode i. The constant 
Wi is related to the gain (^ in the corresponding mode, 

Wi = wg{uji)/eRpg. 

We restrict ourselves to a steady state solution. Elim¬ 
inating the equilibrium population inversion density 
A/’(r), we get the following set of equations for the equi¬ 
librium mode populations rij. 


^-7i -b eRpgW^ 


rdr ’ipfir) 

J V w + J2j WjUjip'jif) 


rii = 0. 


( 20 ) 


Vanishing nds correspond to non-excited modes. Let 
there be N excited modes, We assume 

that we are not far beyond threshold, so that w ^ 
3,nd we may expand the denominator in 
Eq. (| 20 [). We arrive at the following system of linear 
equations (fc = 1 ,... iV) 
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1 


N 


eR 


Po 


1=1 




( 21 ) 


subject to a constraint > 0. Coefficients Ai^,^ are 
given by 

^ jdr'tljf^{r)iPl{f). ( 22 ) 

They are self-averaging quantities with negligibly small 
fluctuations around their mean = l + 25i^i ^, which 

follows from the independent Gaussian distributions for 
[^. Because the correlations between ^ij.i,’s and 
7 i^’s are also negligibly small, we may substitute = 
1 -I- ill Efl- (1^)- Without loss of generality we can 

assume that < ...< -Hj, j giuji^). 

Inverting the matrix we find 


^ _ lik , 1 7»i 

eRj,, N + 2 2gM 2(7V + 2) ^ 5 ( 0 ;,,)’ 

(23) 


The number of excited modes N is restricted by the re¬ 
quirement that all rii^. ’s should be positive. A necessary 
and sufficient condition is 


{2 + N) 


An 

gi^iw) 


E 


Ai 

g{^ii) 


< 2. 


(24) 


Eq. (Q) can be used to determine the probability 
distribution of the number of excited modes, using the 
Porter-Thomas distribution (^) for the statistics of de¬ 
cay rates ji. In the region of parameters where (N) 1 

this mean value can be found analytically from the con¬ 
tinuous approximation of the condition (p4|), 




da a{a) 


^max 



da aa(a) 


2, (25) 


with {N) = da a{a). The density a-{a) of the vari¬ 

ables ai = ^i/giuji) is given by 


cr(a) 



•,x 

du} p{uj)Pi [ag(w)/ 7 (w)] 5 (w)/ 7 (w) 


(26) 


It follows from Eq. ( |^ ) that {N) = Ci{U,/5 liJq)z, where 
z satisfies a cubic equation 

z"^ + \Ci{VI/5u}q)z^ = e. (27) 

To leading order in l/(Af) the term z^ can be neglected, 
which yields a simple answer 

(N) = 3^/3((7iII/<5u;o)^/^e^/^ = (28) 


The general form of this result for any z/ can be derived 
in a similar way, leading to (IVnc) = Cu(^/du}o)e''/'^, 


(N) = (ly + 2 )‘^A>y+‘^) ^ These results are in¬ 

dependent of 7 * as long as 7 * <C 70 - 



FIG. 3. Average number of excited modes (N) versus di¬ 
mensionless pumping rate e (same parameters as in Fig. 1). 
The solid lines are the analytical result (p7|), the data points 
are a Monte Carlo average. The main plot corresponds to 
v = 1, M = 10^ (circles), M — 10® (squares), M — 10"* (di¬ 
amonds). Dashed lines represent the average number (Nnc) 
of non-competing modes. The inset shows the case v = 2 for 
M ^ 1. Note a drastic reduction in the number of excited 
modes. 

To test numerically the analytical results for (N), we 
did a Monte Carlo average over the Porter-Thomas dis¬ 
tribution. For each of 2000 realizations, we ordered the 
modes in increasing order of the ratio loss over gain and 
found maximal N satisfying Eq. (|^). Results for {N{e)) 
are in excellent agreement with the continuous approxi¬ 
mation down to (N) ~ 1 (Fig. 

V. CONCLUSION 


To summarize, we have considered lasing of a chaotic 
cavity coupled to the outside world via v small holes. We 
assumed that the broadening of the cavity modes (due to 
leakage through the holes and absorption by the cavity 
walls) is less than their spacing and used a simple cri¬ 
terion “modal gain > modal loss” as the condition for 
a given mode to be excited. Natural unit of the pump¬ 
ing rate Rp^ is defined such that “maximal gain = mean 
loss.” Because of strong fluctuations of modal widths, the 
probability of lasing can be significantly large for much 
weaker pumping rates than . The distribution of the 
lasing threshold turns out to be wide, with the mean 
much less than Rp^. We have described the multi-mode 
operation as a result of spatial hole burning and found 
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that the average number of excited modes is proportional 
to the power vjiy + 2) of the pumping rate. 
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